Encounter times in spatially-heterogeneous populations and their role in epidemic spread 
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We develop an analytical method to calculate encounter times of two random walkers in one dimension when 
each individual is segregated in its own spatial compartment and shares with its neighbor only a fraction of the 
available space, finding very good agreement with numerically-exact calculations. We model a population of 
susceptible and infected territorial individuals which have such spatial heterogeneity, and which may transmit 
an epidemic when they meet. We apply the results on encounter times to determine analytically the macroscopic 
propagation speed of the epidemic as a function of the microscopic characteristics: the confining geometry, the 
animal diffusion constant, and the infection transmission probability. 
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The spatial propagation of an epidemic in a population is 
highly dependent on the transmission dynamics, in particular, 
the frequency of encounters between individuals, or contact 
rate, and the probability of transmission upon such encounters 
ifTTl . These, in turn, are affected by individual mobility Q, 
by the heterogeneity of the environment [3 ], and by what is 
transmitted, e.g., an infection, a rumor, or information, which 
influences the network of contacts over which the epidemic 
will spread 

Recent studies on the importance of individual-level pro- 
cesses on disease propagation have reemphasized the need 
to develop models that go beyond the assumptions of well- 
mixed, homogeneous populations (SHE), and to link 'macro- 
scopic' features of the spatial spread of an infection to 'mi- 
croscopic' characteristics of the agents carrying the infection. 
In this Letter, we make a considerable advance in that direc- 
tion by focusing on the role that the spatial organization of a 
population plays for the spread of an epidemic through that 
population. We relate the propagation speed of the epidemic 
to pairwise interaction events, consisting of direct encounters 
between one susceptible and one infected individual. 

The rate of encounters between individuals has been ana- 
lyzed in the case where different individuals occupy the same 
spatial region l9lffT1i. based upon recent theoretical develop- 
ments on first-passage times 1 12] in confined geometries lf]~3l- 
[T5ll . Here, we develop a framework to study encounter times 
when the regions occupied by individuals are distinct. 

Model:- Individuals are modeled as random walkers con- 
fined to regions of equal size; a fraction of each region over- 
laps those of the neighboring individuals. This spatial ar- 
rangement is common in territorial animals: individuals are 
confined within a home range fT6lL having exclusive access to 
certain core areas, the territories, 1 17 ], but also sharing other 
regions, the home-range overlap, with their neighbors. 

When a population is arranged spatially, as depicted in 
Fig.[T] on a one-dimensional lattice, infection spread proceeds 
by contact between one infected individual and its neighbor. 
Note that encounters may take place only within the corre- 
sponding overlap region. To represent an SI (susceptible- 
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FIG. 1. (Color online) (a) Three discrete random walkers, with posi- 
tions W\ , W2 and W3 on a one-dimensional lattice, in their respective 
habitats, H\ , H2 and #3 . (b) Mapping of W\ and W2 to the shifted co- 
ordinates (w\,W2) of a single 2D random walker. Static targets cor- 
responding to the locations where the two walkers meet are shown 
for the cases m = 5 (squares), m = 1 (triangles) and m = I (circles). 



infected) model, we suppose that, upon encounter, transmis- 
sion of the infection takes place with probability p. 

The discrete position of the /th individual is denoted by Wf, 
it is restricted to move inside its home range region Hi by re- 
flection at the boundaries, and has overlap regions Ri-\ and 
Ri with, respectively, its left and right neighbor. The sizes of 
both types of region are constant: the Hi each occupy £ sites, 
and the Ri occupy m sites. Interestingly, a similar geometry 
has been independently used in a model for heat conduction 
in one-dimensional systems fT8l[T9) . 

Encounter times:- We first consider the case in which the 
probability p of transmission upon encounter is 1 . To calcu- 
late the mean first-encounter time (MFET) between the two 
walkers, we map their movement onto that of a single random 
walker in 2D, whose allowed positions are represented by the 
vector w = {w\,W2)\ see Fig.[TJb). Here, w; := Wi — (i — 1)A 
is the displacement of walker i from the leftmost site of its 



2 




FIG. 2. (Color online) Unfolding procedure for / = 15: (a) original 
system with reflecting boundaries; and (b) unfolded system. Filled, 
colored sites correspond to targets for the 2D walker for different val- 
ues of m; the cases m = 1, m = 7 and m = 12 are shown by triangles, 
squares, and circles, respectively, (c) Approximation for ra = 12, re- 
placing the region between the inner rhombus and outer square by an 
annulus with inner radius a and outer radius b. 



habitat Hi, with < w,- < £ — 1, and A := £ — m is the distance 
between the centers of Hi and Hi+\, i.e., the mean spacing 
between walkers. 

The possible locations where the two walkers may en- 
counter one other then become a sequence of static target po- 
sitions on a 2D lattice. When m = 1, for example, the overlap 
region consists of only a single site where the walkers can 
meet, and so there is a single target in 2D, namely the bottom 
right corner, w = {£ — 1,0). Asm increases, the set of target 
sites grows, always being arrayed along a diagonal, until the 
limiting case m = £, when the two walkers occupy the same 
space and the targets fill the main diagonal. 

The first-encounter time is now equivalent to the first- 
passage time to the set of targets in the 2D domain, with re- 
flecting boundary conditions at the outer edge of Fig. [IJb). To 
compute this, it is more convenient to consider a symmetrized 
version of the problem: instead of reflecting back a walker 
which attempts to leave the domain, it can be allowed to move 
unimpeded into a reflected copy of the original domain. Doing 
this for all domain borders gives rise to an unfolded version of 
the system ifTTIl : see Fig.[2|a-b). 

When m = 1, the extended system has a "thick" target, i.e., 
four neighboring sites at the center of the lattice, while for 
other values of m a larger rhombus of target sites is obtained. 
The unfolding now allows us to estimate the first-encounter 
distribution of the walkers by computing the first-passage dis- 
tribution to the target rhombus in the continuum limit. 

When the 2D walker starts inside the rhombus, correspond- 
ing to W\ > W2, i.e., with walker W\ to the right of W2 in- 
side the overlap region, one writes the solution of the 2D dif- 
fusion equation for the probability distribution P m (x,y,t) in 
Cartesian coordinates with zero Dirichlet boundary conditions 
on the edge of the rhombus: P m (Q,y,t) = P m (L,y,t) = and 
P m (x,0,t) = P m (x,L,t) — 0, where L = m\fl. The case of 
uniform initial conditions for both w\ and W2 corresponds to 
Pm(x,y,fy = l/^ 2 ' Pm(x, y,t) then has a solution in terms of 
a double infinite series (see, e.g., (TT1 (20|), with each series 
being associated to one of the two orthogonal directions. 

From that expression it is straightforward to calculate the 



survival probability S[ n (t) — f$dx Jq dyP{ n (x,y,t), and finally 
the first-passage probability density F[ n (t) = —dS- m {t)/dt, as 
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The mean first-encounter time is then J^t F{ n (t) dt, giving 
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valid for uniform initial conditions inside the rhombus, where 
the diffusion coefficient D = \, and soC^O.2116. 

When the 2D walker instead starts outside the rhombus, 
we replace the true geometry of the problem in the contin- 
uum limit by an annulus bounded by two concentric circles, 
as shown in Fig.[2jc): one absorbing, of area equal to that of 
the inner rhombus, and hence radius a = with zero 

Dirichlet boundary condition; and one reflecting, with area 
equal to that of the unfolded square, and so radius b = 2£j>jK, 
with no-flux (Neumann) boundary condition. 

We now seek the radially symmetric solution P ou t( r ,t) of 
the diffusion equation in this annulus, where r is the ra- 
dial coordinate, with boundary conditions P ui\r=a = 

and 

d Pout I dr\r=b = 0- F° r reasons that will become clear when 
we discuss epidemic spread, we consider uniform initial con- 
ditions inside the annulus c < r <b, where c > a. Using sep- 
aration of variables in the radially symmetric diffusion equa- 
tion, one can determine the probability of occupation P vx{r,t) 
l20lL the survival probability S O ut(0 = fa r Pontic t)dr, and fi- 
nally, by differentiation, the first-passage probability F ou t(0 = 
—dS 0Ut (t)/dt to reach the inner absorbing circle: 
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with z := f , w := - c and z < w < 1. Here, Ji(s) and Yi(s) 
are Bessel functions of order / of the first and second kind, 
respectively; a n = cc n (z), which depends explicitly on z, is 
the nth positive root of the function <j) z (x) := Jq{xz)Y\{x) — 
Yq{xz)J\{x)\ and p n (s) := Yi(s)Ji(<Xn) - Ji(s)Yi(a n ). Note 
that the roots a n and the functions (f) z and p n all depend on 
the geometry of the annular region, via the parameter z. The 
particular case of uniform initial conditions in the whole annu- 
lar region corresponds to c = a, giving for the corresponding 
first-passage time 
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Combining the two previous results |2| and (|5](, we obtain 
the overall MFET for the 2D walker: 



(^enc) — Pin (^enc)i n ~l~ Pout (^enc) 
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where Pi n (resp. P out ) is the probability that the walker starts 
inside (resp. outside) the target square. When the initial con- 
ditions for the 2D walker are uniform inside the available do- 
main, Pi n is the proportion of sites inside the rhombus, equal 

2 

to ^ m me continuum limit, and P out = 1 — Pi n . We thus have 
a general expression for (T enc ) /^ 2 , which depends only on the 
single dimensionless parameter ™ . 

When z = in the continuum limit, the 2D walker in the 
circular domain is subject only to the reflecting boundary at 
r = b, and the coefficients a n reduce to the roots /5 n of J\ . 
One can then show that for z <C 1 it has the perturbative form 
a n (z) — /} n + Yn(z), where y n (z) decays to zero than slower 
any power of z- For the case n = 1, we have j5\ = 0, giv- 
ing ai(z) ^ V^/He~ 3/4 /z) and Gi(z) ~ [2\n(e^l A /z)}~\ 
whereas for n > 1, we have % > and G n (z) ~ ln _2 (z), 
with the ln(z) dependence of G n (z) resulting from the small- 
argument expansion of the Bessel function of the second kind 
of order zero. For small z, the n = 1 term in the series thus 
increases logarithmically, whereas the terms with n > 1 decay 
to zero. 

We thus obtain the following approximation for the mean 
encounter time: 

(Tenc) Cm* 8, ( e-y 4 V2( \ 

which exhibits the well-known logarithmic dependence of 
mean first-passage times for small target size [Q21 US . Note 
that from the exact first-passage distribution §5§ to a circular 
target, we can determine all higher-order moments of the first- 
encounter time. These results, together with a detailed anal- 
ysis of the dependence on the ratio between the radius of the 
confining domain and the target, will be reported elsewhere. 

We test the accuracy of these results by comparing them 
with numerically-exact values, obtained by solving sparse 
systems of linear equations, consisting of recurrence relations 
for the mean first-passage time to reach a target from a given 
site in terms of the same quantity starting from the neighbors 
of that site. Figure [3] shows numerically-exact results for sev- 
eral combinations of m and £, which all collapse onto a single 
curve; this, in turn, is very well approximated by the analytic 
expression ([6]) for the continuum limit. We also compare the 
first-order approximation ([7]), which mirrors the analytical re- 
sults for small j, but is not accurate for ™ > 0.3. 

Epidemic spread:- To pass from the microscopic interac- 
tions between neighboring walkers to the macroscopic propa- 
gation of an epidemic, we consider a chain of habitats Hi, each 
with one confined walker, i = 0, 1, . . . , as shown for three 
walkers in Fig. [T] Starting with the left walker infected, the 
time to infect walker n to the right is given by = YH=\ x u 
where T ; is the time for walker i — 1 to infect walker i, which 
are independent random variables. 

Each pair of walkers i — 1 and i may be treated as an identi- 
cal copy of the same system (except the leftmost pair). How- 
ever, the initial conditions are no longer uniform over the en- 
tire square, since when Wi is infected by its left neighbor W/_i, 
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FIG. 3. (Color online) Dependence of (T enc ) /£ 2 as a function ofm/£. 
Numerical results are shown with (black) dashed lines, the analyti- 
cal continuum limits {6} with a (red) solid line, and the approximate 
expression {7} with a (blue) dash-dotted line. The inset shows the 
corresponding curves on a log-log plot. The mean first-encounter 
time relevant to epidemic propagation (walker W\ uniform in its left 
overlap region, i.e., restricted initial conditions) is shown only for 
j < \, delimited by a vertical dotted line, since we consider only 
cases for which individuals always maintain an exclusive area. 



its initial position must lie inside the overlap region R[-\ with 
that neighbor. For the 2D walker, this corresponds to taking 
a distribution that is uniform only within two vertical strips, 
of width m and height 2i, at the left and right edges of the 
unfolded square domain shown in Fig. |2jb). Since, by sym- 
metry, the same result is also obtained starting from two hor- 
izontal strips, this can be approximated by a square strip of 
width m around all edges of the square domain, neglecting 
double counting in the corners where the horizontal and ver- 
tical strips overlap. We then approximate this by the annulus 
c <r <b, where c = 21(1 — j)/ y/n, so that its area is equal 
to that of the square strip. 

Taking uniform initial conditions in this annulus, and us- 
ing Pi n = in Eq. ([6J, gives the corresponding mean first- 
encounter time as (T enc ) = (T enc ) out , whose expression is as 
in Eq. ([5]), but with G„(z, 1) replaced by G n (z,w), with w — 
m/[V2£(l —m/£)]. This analytical result is also shown, and 
compared to numerically-exact values, in Fig. |3j again with 
excellent agreement. 

If the infection is not transmitted at every encounter, but 
only with a certain probability p < 1, then the two walkers 
may meet repeatedly before the infection is transmitted, at any 
position in the overlap region. We thus must take into account 
returns to the target set, in which the 2D walker starts from one 
of the possible targets and comes back to hit any of them. To 
do so, note that the infection is transmitted at the &th encounter 
with probability pk := p(l — p) k ~ l , so that the mean infection 
time is given by (t) = J^° =1 p^ (z\k) , where (z\k) is the con- 
ditional mean infection time, conditioned on infection occur- 
ring at the &th encounter. Since (z\k) = (T enc ) + (k — 1) (T re t), 
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FIG. 4. (Color online) Position co(t) of the epidemic front, as a 
function of time t, in a ID population of territorial individuals with 
habitat size £ = 20 and overlap size m = 5. Analytical predictions 
(dashed lines) and numerical results (thick lines, averaged over 100 
realizations) for the mean position are shown for different transmis- 
sion probabilities p. The time scale shows the number of simulation 
steps; (T e nc) — 629 for this choice of £ and m. Trajectories for several 
realizations are shown in light lines: (red) solid for p = 1; (green) 
dotted for p = 0.1; and (blue) dashed-dotted for p = 0.01. Small 
crosses indicate the analytical predicted spread of possible front po- 
sitions at a given time for p = 1 . 

where (T ret ) is the mean return time, we obtain 

(T) = (T enc ) + (i-l)(T ret ). (8) 

The mean return time (r ret ) to a set U can be calculated for 
discrete stochastic models, using the Kac recurrence lemma 
EniED, as (T ret ) = |^|/|^|, where £1 is the set of all possible 
configurations and | • | denotes the size (number of configura- 
tions) of a set. Taking U as the set of targets, we have \U\ = m, 
giving (T ret ) /£ 2 = 1/m, which is no longer a function of the 
dimensionless parameter mji. Using this in Eq. ([8]), we obtain 
the mean infection time (t) for any p < 1. 

We can now put the previous results together to calculate 
the propagation of an epidemic in the chain of habitats. We are 
interested in the position co(t) of the epidemic front, i.e., the 
position of the infected walker which is furthest to the right 
l22l . The nth walker is infected at time T^, at which time 
the front is at co(tW) = n(£ — m). Taking the average over 
realizations, the propagation of the front is thus linear in n, 
with mean speed (£ — m)/(z). 

Figure [4] compares the analytical result for the mean front 
position with numerical simulations, averaged over 100 re- 
alizations, showing excellent agreement for different values 
of the transmission probability p\ several individual trajec- 
tories are also shown. For p = 1, we additionally plot the 
analytically-predicted range of front positions at a given time, 
using the so-called law of the iterated logarithm 1231 , which 
shows that the infection time of walker n should be within the 



interval n (t) ± y/2n loglogftcr(T) for large enough n. Here, 
<j(t) denotes the standard deviation of the first-encounter 
time, which we have calculated analytically from the exact 
expression ^ for F out . Again this agrees with the numeri- 
cal results, shown as a cloud of realizations of the epidemic, 
which are contained inside these bounds for large n. 

In summary, we have studied a model of suscepti- 
ble/infected territorial animals in ID with overlapping home 
ranges. We have shown how to calculate first-encounter times 
between two neighboring individuals, finding excellent agree- 
ment between analytical and numerically-exact results. Tak- 
ing into account that the infection may not be transmitted at 
each encounter, we have also obtained analytical estimates for 
the mean infection times. We have successfully applied these 
results to infer the macroscopic propagation speed of an SI 
epidemic, as a function of the heterogeneous spatial distribu- 
tion and the interaction dynamics between individuals. 
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